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1. Introduction 



One of the foremost candidates for dark matter are WIMPS, particles with masses (hundreds 
of GeV) and interactions in the weak force range. Particles with these properties are generated 
thermally in the early universe but decouple from the other particles early enough such 
that a sizable number of them survived at present times. If they are stable their mass and 
density would have the correct order of magnitude to be the dark matter observed in rotation 
curves of galaxies, galactic clusters, and cold dark matter models of structure formation. 
Recently, observations of the cosmic ray flux showed an excess in the high energy component 
as compared to the standard models [1] describing the diffusion of cosmic rays originating 
from supernova remnants [2, 3, 4, 5, 6]. Furthermore, the enhancement is observed in the 
electron flux but not in the proton flux. One possible explanation to the anomalous cosmic 
ray data is that it is the result of dark matter particles annihilating [7, 8, 9] into light 
particles. If these light particles have a mass below hadronic scales, they will decay mainly 
into electrons which would explain why the enhancement is not seen in the proton flux. One 
problem with this interpretation of the data is that, given the known dark matter density 
and the weak scale annihilation cross sections, the electron flux falls short from the observed 
flux by one to three orders of magnitude (in the absence of significant dark matter clumping 
and dependent on the observational data set considered). It has been advocated [10] that 
models including a long range attractive force between dark matter particles can enhance the 
annihilation cross sections to the observed levels, a phenomenon refered to as the "Sommerfeld 
enhancement" [11, 12, 13, 14, 15]. If this force arises from the exchange of a particle with a 
mass smaller than hadronic scales but larger than the electron mass scale, the standard success 
of Big-Bang nucleosynthesis is preserved. This light particle (referred here on as </>) can also 
be the product of dark matter annihilation, nicely explaining the enhancement of the electron 
but not the proton flux. If this picture is correct, the cosmic ray data should have a peak 
around the mass of the dark matter particle. This is seen in the older data but not in the 
most recent GLAST results [6]. An enhancement caused by the exchange of a light boson 
is also the essential ingredient of a unified theory for dark matter[10, 16] that also addresses 
other observations like the WMAP haze and the DAMA results. 

Assuming that the picture above is correct, one has to deal with the fine tuning required 
for the presence of a particle with mass in the GeV range in a theory with a typical scale 
in the TeV range. Although solutions to this problem involving supersymmetry exist (for 
example, see Ref. [17, 18, 19]), a simple, natural solution would be to assume that the 4> 
particle is a pseudo-Goldstone boson. The small mass of <j) would be, in this case, "technically 
natural" in the sense that radiative corrections would be zero in the = limit. The force 
generated through the exchange of a pseudo-scalar pseudo-Goldstone boson is, in general, 
very complicated and spin dependent. In fact, it is a close analogue to the force between 
two nucleons generated by pion exchange. Let us assume the dark matter particles are spin 
1/2 Dirac fermions. The force due to 4> exchange, just like the nuclear force due to pions, 
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contains two pieces. One is a central force of the Yukawa type 1 . The other, the tensor force, 
operates only on spin triplet states and mixes L = and L = 2 components of the wave 
function. At short distances, r <C l/m^, the tensor potential goes as ~ 1/r 3 and it becomes 
the dominating piece of the interaction. An observer might naively think that a 1/r Yukawa 
potential will be dominant over a larger range than a 1 /r 3 tensor potential. However, this is 
not the case. When comparing the form of the Yukawa potential to the tensor potential, 
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one finds that the Yukawa could only dominate when r > However, both potentials 

are exponentially suppressed when r ~ l/m^. Thus, within the range of these potentials, the 
tensor part dominates the Yukawa force. 

An additional complication is that the singular tensor force makes the Schrodinger equa- 
tion ill-defined. Therefore, we are led then to a renormalization procedure to make the 
non-relativistic effective theory considered here meaningful. In the present case, the renor- 
malization program consists of introducing, in addition to the long range force, a short dis- 
tance interaction generated by a contact, four-fermion operator in the low-energy effective 
theory. The range of this interaction is of the order of the inverse mass of other states in the 
theory that have been integrated out, that is, of the order of TeV -1 . The strength of this 
potential can be determined only through a matching calculation to the full theory in the 
dark matter sector and, as such, it is an additional parameter from the point of view of the 
effective theory. The complications due to the mixing of the partial waves and the need for 
the renormalization procedure form the bulk of the present paper. In the next section, the 
issues arising in the calculation of the Sommerfeld enhancement in the presence of a singu- 
lar potential are discussed in the simpler context of a central, single-channel 1/r 3 potential. 
This example is not entirely without physical interest as models of dark matter with electric 
dipole moments would contain such potential. In section 3 we present the results of the full 
coupled-channel calculation. The conclusions are summarized in section 4. 



2. Sommerfeld enhancement with single-channel singular potential 

The Sommerfeld enhancement is a simple quantum-mechanical effect that, at least in the 
simpler case of uncoupled channels, has been discussed in the literature several times. Here 
we briefly review the single channel calculation in order to highlight a couple of points that 
are relevant to the non-central, singular potential case of interest. 

Let us consider the annihilation of a pair of Dirac particle/antiparticle of mass M, attract- 
ing each other with a central, spin-independent potential proportional to 1 /r 3 . This potential 
could be generated by the exchange of an electromagnetic or "dark" photon between particles 

1 For more details on enhancement from a dominant central Yukawa interaction, see Ref. [20, 21]. 
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with a permanent electric dipole. 2 . In addition, there will be a short distance interaction at 
a scale R comparable to the one set by the mass (R ~ 1/M). The inclusion of this short dis- 
tance piece is dictated by the renormalization properties of the singular potential. If we take 
the quantum mechanical potential model as an effective theory of a microscopic relativistic 
theory, this short distance potential arises from a four-fermion operator needed to match the 
high energy theory. 

We look for a solution to the Schroedinger equation satisfying the usual asymptotic 
condition at large separations 
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where r is the separation between the two particles, k the incoming momentum and / = 
J^ (2Z + l)fiPi(cos9) is the scattering amplitude. The wave function can be expanded in 
partial waves as 
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VKf) = £(2Z + l)^l Pl ( cos e), (2.2) 

1=0 

and for central, spherically symmetrical potentials, the Schroedinger equation decouples into 
separate equations for each partial wave: 

-itj'(r) + MV{r)ui{r) = k 2 ui{r). (2.3) 

The solutions to eq. (2.3) behave at large r as free waves 

ui(r) -» Asm(kr -lir/2 + 5i). (2.4) 

The constants A and Si are fixed by the two boundary conditions. One is ui(0) = 0, needed in 
order for ?/>(0) to be finite. The second one comes from the asymptotic condition in eq. (2.1), 
which implies 

JpiSi 

A = -JT- (2 - 5) 

The scattering amplitude is related to the phase shifts <5; through e 2i<5 ' = 1 + 2ikfi and does 
not depend on the normalization factor A. The same is true about the annihilation cross 
section enhancement. This enhancement is given by the ratio between the interacting and 
the free wave functions at origin (assuming a zero range bare annihilation amplitude) 

2 
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Thus, in the single channel case, the normalization in eq. (2.5) plays no role and can be safely 
taken to be A = 1 (both for the free and interacting wave u). In the expression for 

2 Recent studies [22, 23] found bounds on the electric and magnetic dipole moments from experiments. In 
the context of our study, these correspond to (10 -10 a e ) < G < (1 a e ) where a e = 1/137. Bounds on a 
dominant dipole interaction in the dark sector from a massive dark photon have yet to be explored. 
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Figure 1: Annihilation process enhanced by multiple Goldstone boson exchanges. Solid lines represent 
the dark matter particles, dashed lines the Goldstone boson and wavy lines may or may not be the 
Goldstone particle. 

the Sommerfeld enhancement above, it was assumed that the annihilation proceeds through 
a zero range interaction. More realistically, the mechanism of annihilation occurs within a 
distance scale ~ 1/M. This would be true, for instance, if the annihilation goes through the 
process depicted in Fig. (2) as the vertical internal fermion lines are off-shell by an amount 
~ M. The annihilation rate would then be determined by a weighted average of the wave 
function over a region of size ~ 1/M around the origin. However, for well behaved potentials, 
the wave function varies little over this region and eq. (2.6) suffices. 

As mentioned above, there is an additional issue when the potential is more singular at 
origin than 1/r 2 . The Schrodinger equation with this potential is meaningless and leads to the 
"fall-to-the center" scenario due to the extreme short distance attraction. Bound states exist 
with arbitrarily high (negative) energy and the hamiltonian is not bounded from below. The 
way out of this problem is well known. The quantum mechanical model we are considering 
should be viewed as the low energy effective theory of the underlying relativistic field theory 
describing the system. The validity of the non- relativistic description and the 1/r 3 potential 
is limited to distances larger than a certain scale R where the effective theory breaks down 
and the potential for shorter distances cannot be trusted. This fact, of course, does not 
invalidate the use of the Schrodinger equation to describe low energy annihilation. It only 
means that a renormalization procedure should be carried out. In the present case we carry 
out the renormalization program by splitting the effective potential into two pieces. The first, 
valid for r > R, is the 1/r 3 potential generated by a light particle exchange. The second, 
valid for r < R is a constant potential generated by a four-fermion term in the effective theory 
lagrangian. The distance R is not a parameter of the model and identical descriptions of the 
low energy physics should result from different values of R, as long as the low energy couplings 
of the effective theory are made to depend on R appropriately. For non-singular potentials, 
the dimension 6 four-fermion operator describing the short distance interaction is irrelevant 
and appear only at higher orders of a low energy expansion. For potentials more singular than 
1/r 2 this is not true and the short distance potential is necessary for a consistent calculation 3 . 
The height Vo and the range R of the potential are arbitrary. However, the requirement that 
the low energy observables are independent of the regularization procedure (the values of R 

3 Additionally, this short distance potential can have an 0(k 2 ) correction, but for small k, this correction 
is negligible compared to Vo. 
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and Vo) fixes the dependence of Vq on R. The whole arbitrariness of the procedure is then 
reduced to the arbitrariness of one parameter which can be taken to be the value of Vo at 
one specific value of R. Physical results will depend on the value of this parameter. Given a 
microscopic model one can calculate the scattering amplitudes by matching the microscopic 
model to our effective theory. As a consequence, the Sommerfeld enhancement will not have 
the universal character it has when only non-singular potentials appear. 

The fact that the all low energy observables are rendered cutoff independent by adjust- 
ing one short distance coefficient is shown, in the context of singular potential in quantum 
mechanics, in Ref. [29, 24]. In our case, we solve eq. (2.3) with the potential 

We will consider here a purely s-wave bare annihilation amplitude of short range. "Short 
range", in this context, means a range comparable to the cutoff R. In principle it is possible 
to construct models where the range of the short distance elastic interaction, annihilation 
amplitude, and mass are given by different scales. We will only consider here the more 
natural case where both ranges are comparable to R ~ 1 /M. The annihilation amplitude is 
subjected to the same renormalization procedure described for the elastic interaction. The 
bare annihilation amplitude is a cutoff dependent quantity and only the amplitude dressed 
by the potential is physical and cutoff independent. Again, given a microscopic model one 
can evaluate the value of the annihilation amplitude in that model. For our purposes, we 
regulate in the annihilation amplitude in a manner similar to the elastic potential 
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where uo(r) is the wave function in the absence of the potential with the same normalization 
as u(r) at large r. Different ways of smearing the ratio tt(0) /uq{0) lead to very similar results. 
The numerical solution of eq. (2.3) does not bring additional difficulty since the form of the 
wavefunction is fixed by the u(Q) = boundary condition. 

In Fig. (2) we present results for a sample of parameter values. These parameters were 
chosen to mimic the physics of an actual exchange of Goldstone bosons in a theory for dark 
matter as mentioned in the introduction. Thus, the strength of the potential (G ~ 1) is of 
the order of the one generated by Goldstone bosons with a decay constant of order ~ M, 
and the hierarchy between the masses of the dark matter particle and the Goldstone boson 
is TeV/GeV ~ 1000. A general trend towards larger enhancements for larger values of G 
is observed. This is expected as a larger value of G means a more attractive potential. In 
addition there are wild oscillations superposed on this general trend as either G or Vo is 
varied. The peaks in the enhancement occur when there is a bound state close to threshold. 
Let us compare these results to the enhancements generated by other central potentials. In 
the case of the Yukawa potential, analytical estimates are available and show a general trend 
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Figure 2: Sommerfeld enhancement 5 as a function of the potential strength G for two value of 
V = l/MR 2 , W/MR 2 . The other parameter used are M = 1000 GeV, m = 1 GeV, k = 0.2 m and 
R = l/M. 

of S growing proportionally to the strength of the potential. Numerical calculations show 
that, in addition to this trend, there are regions of parameter space where the enhancement 
is much larger. The enhancement due to a square well potential is analytically calculable 
and again shows the large enhancements as bound states cross the threshold. When both 
long and short distance potentials are combined it is difficult to obtain any reliable analytical 
estimates. There is also a complicated dependence on the momentum k. Since the 1/r 3 
potential does not make sense by itself, we always have to consider the short potential with 
it, which makes analytical estimates difficult. 

3. Sommerfeld enhancement though Goldstone boson exchange: singular 
couple-channels 

Goldstone-like pseudo-scalars with derivative couplings to fermions in the dark sector give 
a naturally light exchange particle for dark matter interactions. It would be convenient 
for building dark matter models if such as mechanism could lead to an enhancement of 
annihilation cross sections. This possibility was considered in Ref. [10] but quickly discarded 
since the part of the potential surviving in the m = limit (the tensor force) vanishes when 
averaged over the s-wave initial and final state. While true that the tensor contribution leads 
to a change in the orbital angular momentum of AL = 2, for a fixed total angular momentum 
J, this contribution can take a state with L = £toL = £ + 2 and then back to L = I. Since 
this process can occur an arbitrary number of times, a coupled-channel analysis is necessary in 
order to account for this effect as in the one pion exchange potential in the deuteron channel. 

The exchange of pseudo- Goldstone bosons between spin- 1/2 fermions generates a singular 
potential similar to the one considered in the previous section. In fact, this problem is very 
similar to the scattering of two low-energy nucleons and we refer to the literature [24, 25, 26, 
27, 28] for a more detailed discussion of the power counting for this non-relativistic effective 
theory. In that scenario, the long distance potential is also dominated by the exchange 



of pseudo-Goldstone particles (pions). There are, however, a few differences between the 
present problem and the nuclear case. First, the sign of the one-Goldstone exchange potential 
between particle and anti-particle is the opposite to the one between two particles. Second, 
the separation between the M and m scales is more pronounced than in the nuclear case. 

At the lowest order in the low-momentum expansion the pseudo-scalar Goldstone boson 
couples to the axial current of the dark matter particle, 

f 2 X^'x B -f, (3.1) 

where x 1S the dark matter particle, assumed to be a spin 1/2 Dirac fermion, <p the Goldstone 
boson, g the axial coupling of these fermions and / the Goldstone boson decay constant. 
Naturality arguments suggest a value of g ~ 1 and 4irf ~ TeV, as well as M ~ TeV. The 
constraints mentioned above and discussed in [10] imply in m > 0.3 GeV. Let us examine 
the one-Goldstone exchange between a fermion and anti-fermion in our dark sector with total 
angular momentum J = 1 . The resulting one-pion exchange potential is well known in nuclear 
physics and is given by 

V(r) = Voir) a\- (72 + Vr(r) (3r ■ a\r ■ a 2 - o\- a 2 ) (3.2) 

where 
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and a = g 2 /{IQnf). The central part of the potential corresponds to a AL = transition, 
while the tensor part of the potential corresponds to a AL = 2 transition (the spin structure 
yields a traceless symmetric tensor). In the single channel case, we were interested mainly 
in the s-channel annihilations since the annihilation rates are larger there due to the absence 
of a centripetal barrier. Here, we want to consider the same s-channel annihilation. The 
J = L = S = channel is of no interest since the Yukawa potential in this state is repulsive. 
We are led to consider then the J = 1,S = 1,L = 0,2 coupled channels. The most general 
wavefunction with S = 1 , mj = quantum numbers is 4 



C^%_ l{ r)\l, 1) + Ct^m, "I) + ^ +1)0 V, (r)|l, 0)" 

1=0 ' 

OO / \ 

(3.3) 

4 There is in addition a J = I component that, however, does not couple to the J = £—l,£-\-l components 
in the absence of a spin-orbit force. 
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where the Clebsch-Gordan coefficients, C^ J , are given by 



n (e+i)oe _ I I r (t+i)oi _ I 1 n (i-i)oe _ / 1 + 1 

U±1 ~ V 2(2^ + 1)' °° ~V2£ + 1' ±x "V 2(2*+ 1)' G ° " V2^+l' 

(3.4) 

and the kets represent spin states \S, mg). The asymptotic boundary condition analogous to 
Eq. (2.1) when the initial particles are in a mj = state is 

^(r)^e^|l,O} + /(0,<7)— = ^f;(2€+l)P,(cos0)[ e ^-(-l)V^]|l,O) + /(f )( 7)^ 

£=0 

(3.5) 

where <r represents the spin degrees of freedom. For distances larger than the potential range 
the wave function is free and approaches 

( „— iei ( ikr+2t€( / -i \£ —ikr 



u e (r) -> sin(A;r - Ztt/2 + e*) = — '-e~ lte e lkr+Met - {-ife 

2i V 

w t (r) - Q sm(Ar - Ztt/2 + 0) = cj^-e^ (e ikr+2i <* - (-1) V*^) (3.6) 

Matching the incoming spherical wave part of Eq. (3.5) and Eq. (3.6) we arrive at the relations 

jt-piu i^e^ 1 r- 

B e = —— Vtoy/i + l, C e = — \/W£ (3.7) 

k k 

Like in the single channel case the overall normalization of the wave function is unphysical and 
cancels out in the computation of the Sommerfeld enhancement. The relative normalization 
of the ug and we components, however, is meaningful. To make things worse, it appears as 
if the knowledge of the phase shifts ei and ^ is needed before the Scrodinger equation can 
be solved. Fortunately, there is a simple algorithm to construct the wave function with the 
proper boundary conditions. Let us restrict ourselves to the J = l,mj = case of interest. 
Then only no and w 2 mix under the influence of the tensor force (from now on we drop the 
subscript on u and w). We first find two linearly independent solutions (u,w) and (u,w) 
satisfying it(0) = w{0) = and u(0) = w(0) = 0. At large r, they will behave like free waves 
with arbitrary normalizations: 

u(r) -> Aie~ ikr - B x e ikT u(r) -> A x e~ ikr - B x e ikr 
w(r) -> A 2 e~ ikr - B 2 e ikT ' w(r) -> A 2 e~ ikr - B 2 e ikr ' 

It is easily verified that the relevant linear combination 

(u(r) u(r)\ (a 1 AiY 1 ( 1 \ 
" \w(r) w(r)J \A 2 A 2 J \-^2) { - ] 

satisfies the desired boundary condition eq. (3.5). From the solution in eq. (3.8) we can 
calculate the Sommerfeld enhancement. Assuming that the short distance operator describing 
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Figure 3: Logarithms of the Sommerfeld enhancement S as a function G for Vq 
R = 1/M and coupling is a = G/M 2 , where M = 1000 GeV, and m = k = 1 GeV. 

the annihilation is dominated by the s-wave contribution we calculate 

f R R dr u 2 {r)e- r2 / R2 



-W/(MR 2 ), 



s 



(3.9) 



The resulting Schrodinger equation mixes only states with the same J. In the case of the 
J = 1 state we have 



where 



-*"(r) + M(V L (r)9(r - R ) + V s (r)9(R - r))*(r) = k 2 ^(r), 
u(r)\ T . . . / Vc(r) ^81/ T (r) N / 



V8V T (r) Vc(r) - 2V T (r) + 6/r 2 







(3.10) 



. (3.11) 



Notice that the short distance potential arises, in the effective theory, from a four-fermion 
interaction with no derivatives contributing only to the s-wave interaction. Short distance 
d-wave contributions exist but are of higher order in the low energy expansion. It is also 
worthwhile to note that the reversed sign of the potential, as compared to the nuclear case, 
does not imply that the Sommerfeld enhancement will not occur. In fact, the potential matrix 
above when diagonalized yields one repulsive and one attractive eigenvalue. Upon inspection, 
at small r, the repulsive eigenvalue is suppressed compared to the attractive one. 

Examples of numerical results for the Sommerfeld enhancement are shown in Fig. (3) 
as a function of the potential strength G and in Fig. (4) as a function of the short distance 
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Figure 4: Sommerfeld enhancement S as a function of the short distance potential Vq. For the 
coupling a = l/M 2 , R = 1/M, where M = 1000 GeV, and m = k = 1 GeV. 



potential Vq. They exhibit a similar oscillatory pattern as the single channel calculation in 
Fig. (2). The enhancement is large only at certain combination of parameters where a bound 
state close to threshold exists. Contrary to the single channel 1/r 3 potential, there does not 
appear to exist any trend towards larger enhancement for stronger potentials (larger G). In 
this sense the coupled channel problem behaves like a short distance potential (keeping in 
mind that the oscillatory behavior of the enhancement depends on G) . 

The periodicity of the enhancement as a function of \/— MVqR shown in Fig. (4) is easy 
to understand. The wave function for r > R depends on the value of i/j(R) and ip'(R) only, 
not on the behavior of ip{r) at distances smaller than the cutoff R. The wave function inside 
the square well, u{r) ~ sin(rV — MVq + A; 2 ), exhibits this periodic behavior for small values 
of k 2 . As a consequence, the value of R needed to keep the low energy physics fixed varies 
periodically with Vq, a feature well known in non-relativistic effective theories [29]. While 
most Vq values lead to moderate enhancement (S ~ 1 — 10), there exists a significant range 
of Vq values that can lead to S ~ 10 — 1000 and several near resonance points that can lead 
to S > 1000. Since the quantity Vq is dependent on the higher energy physics that has been 
integrated out, the magnitude of these enhancements via pseudo-Goldstone exchange is not 
a universal quantity and can only be determined within the context of a specific microscopic 
model through the matching of the low energy effective theory to its ultraviolet completion. 
The outcome of this matching is the value of Vo (given a value of the regulator parameter R) 
that can then be used to determine the enhancement S as explained above. 
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4. Discussion 



We discussed the enhancement of the annihilation rates of dark matter induced by the ex- 
change of a pseudo-Goldstone boson. From the point of view of model building for a dark 
sector, this choice has the advantage of naturally generating the required hierarchy between 
the WIMP mass and the mass of the particle mediating their interaction. For the case of 
spin-1/2 Dirac dark matter that we considered, we find modest enhancement for most of the 
parameter space. Still, only moderate fine tuning is necessary to achieve the enhancement of 
order 10 suggested by the more recent FERMI data. The enhancements produced are qual- 
itatively similar to the enhancements generated from a short-distance potential, and hence 
can be large only with a certain amount of fine tuning . It is a thorny issue to quantify how 
much fine-tuning is required for a certain enhancement. It is not a priori clear how to assign 
probabilities in the effective theory parameter space. For instance, should the probabilities 
be uniform in Vq or \/Vo ? 




We have not considered any details of model building of the microscopic theory of dark 
matter that could realize the Goldstone exchange enhancement we discussed. We point out 
though that an obvious model with the necessary features would be a scaled up version of 
QCD in the dark sector. The analogy with nuclear physics are then very strong. WIMP's 
would be the "baryons" and the (ft particle would play the role of the "pions" . Assuming the 
existence of more than one light "quark" , we expect a splitting between different "baryonic" 
states to be of the order AM ~ in 2 /A (A ~ M ~ 4irf is the scale of this theory). For 
Awl TeV, m & 250 MeV , AM « 0.1 - 1 MeV range. WIMP excited states with similar 
splittings have been invoked in the "Exciting dark Matter" [30] scenario also incorporated 
in the unified description advocated in [10]. A similar QCD-like theory was proposed in 
Ref. [31], but, unlike our model, baryon formation is suppressed leaving the dark mesons as 
the primary dark matter candidate. 

Another question not investigated here but that deserves attention is the possible en- 
hancement of the annihilation rates through the formation of bound states over cosmological 
times previous to their decay into light standard model particles [32, 33]. In the absence of 
very large or infinite range repulsive forces (the analogue of the Coulomb repulsion in nu- 
clei), the formation of very large "dark nuclei" seem to be possible. This lumpiness of the 
distribution of dark matter will have a strong effect on any estimate of their annihilation 
rates. 
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